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ABSTRACT 

The perturbation caused by planet-moon binarity on the time-of-arrival signal 
of a pulsar with an orbiting planet is derived for the case in which the orbits of 
the moon and the planet-moon barycenter are both circular and coplanar. The 
signal consists of two sinusoids with frequency (2n p — 3nb) and (2n p — rib), where 
n p and n b are the mean motions of the planet and moon around their barycenter, 
and the planet-moon system around the host, respectively. The amplitude of the 
signal is equal to the fraction sin I[9(M p M m ) /16(M P + M m ) 2 }[r / R} 5 of the system 
crossing time R/c, where M p and M m are the the masses of the planet and moon, 
r is their orbital separation, R is the distance between the host pulsar and planet- 
moon barycenter, / is the inclination of the orbital plane of the planet, and c is 
the speed of light. The analysis is applied to the case of PSR B1620-26 b, a pulsar 
planet, to constrain the orbital separation and mass of any possible moons. We 
find that a stable moon orbiting this pulsar planet could be detected, if the moon 
had a separation of about one fiftieth of that of the orbit of the planet around 
the pulsar, and a mass ratio to the planet of ~ 5% or larger. 
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Subject headings: planetary systems - - pulsars: general - - pulsars: individ- 
ual(PSR B1620-26) — stars: oscillations 



Introduction to Extra-solar Moons 



In the past decade and a half, over two hundred and fifty extra-solar planets have 
been discovered^ With the data expected to be produced by satellites such as COROT 
( Auvergne et al.||2003 ) and Kepler ( Basri, Borucki fc KockpOOB" ), it will not only be possible 



to find smaller planets, but moons of those planets as well (Szabo et al. 2006). As a result 



the detectability of extra-solar moons is starting to be explored in terms of their effect on 



planetary microlensing (Han & Han 2002) and transit lightcurves (Sartoretti & Schneider 



1999 Szabo et al. 2006). Upper limits have already been placed on the mass and radius of 



putative moons of the planets HD 209458 b (Brown et al. 2001), OGLE-TR-113 b (Gillon 



et al. 2006) and HD 189733 b (Pont et al. 2007). 



While the limitations of microlensing and the transit technique for detecting moons 
have been discussed and used in the literature, the limitations of other techniques such 
as the time-of-arrival (TO A) technique have not. This technique involves determining the 
variations in line-of-sight position to the host star, usually a pulsar, using the observed time 
of periodic events associated with that host. The aim of this analysis is to explore exactly 
what the TOA signal of a planet-moon pair is, and relate it to the planetary systems that 
can give the most precise timing information, the systems around millisecond pulsars. 



2. Review of Planetary Detection Around Millisecond Pulsars 

The first planetary system outside the Solar System was detected around the millisecond 



pulsar PSR 1257+12 (Wolszczan & Frail 1992). This detection was made by investigating 



periodic variations in the time of arrival of its radio pulses using a timing model. An example 
timing model for the case in which the planet's orbit around the pulsar is circular is: 



(t N - t ) = (T N - T ) + AT corr + TOA per%p (M s , M p , R, I, 6 (O)), 



(1) 



where to and are the times at which the initial and N th pulses are emitted in the pulsar's 
frame, T and T/v are the times the initial and N th pulses are received in the observatory's 
frame, and the term AT corr acts to change the frame of reference from the observatory on 



1 See, for example, 



http://cxoplanet.eu/catalogue.php 
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Fig. 1. — Coordinate system used in the analysis of the TOA perturbation caused by a moon. 
The subscripts s, b, p and m refer to the pulsar, planet-moon barycenter, planet and moon, 
respectively. The diagram on the left shows the quantities used to describe the position of 
the three bodies in their mutual orbital plane, while the diagram on the right shows the 
relationship between this orbital plane and the observer. 



Earth to the barycenter of the pulsar system (see |Backer|1993| for a more complete discussion 
of the components of AT corr ). The final term represents the effect of a planet on the motion 
of the pulsar, where R is the planet-pulsar distance, I is the angle between the normal of 
the planet-pulsar orbit and the line-of-sight and </>&(0) is the initial angular position of the 
planet measured from the x-axis, about the system barycenter. In addition, M s and M p are 
the mass of the pulsar and the planet, respectively. 

Currently, four planets around two millisecond pulsars have been discovered, three 
around PSR 1257+12 QWolszczan fc Frail|1992||Wolszczan|1994| and one around PSR B1620- 



26 (Backer, Foster & Sallmen 1993). These four planets include one with mass 0.02 Earth 



masses, the lowest mass extra-solar planet known. This low-mass detection threshold, in 
addition to measurements of orbital perturbations (e.g., the 2:3 orbital resonance between 
PSR 1257+12c and PSR 1257+12d, |Konacki fc Wolszczan|2003| , demonstrate the sensitivity 
of the TOA technique. Indeed, as a result of their high rotation rate (and thus large number 
of sampled pulses) and low level of noise activity, millisecond pulsars make optimal targets 



for high precision TOA work (see e.g., Cordes 1993). 
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3. What is the TOA Perturbation Caused by a Moon? 

In order to investigate the perturbation caused by planet-moon binarity, the timing 
model presented in equation ([!]) must be updated to include effects due to the presence of 
the moon. For simplicity, we consider here only systems in which both the orbit of the planet 
and moon around their common barycenter, and the orbit of the planet-moon barycenter 
around the pulsar, are both circular and lie in the same plane. An example model taking 
this assumption into account is: 



N 



N 



T ) + AT corr + TOA pertiP (M s , M p + M m , R, 1, 0^(0)) 

+ TOA pert , pm (M s ,M p ,M m ,r,R,I,M0)A P (0)). (2) 



We have explicitly modified TOA peTt p to indicate that it depends on the combined planet- 
moon mass, and included another term, TOA perttPm , to account for planet-moon binarity. 
Here M m is the mass of the moon, r is the distance between the planet and the moon, and 
4> p (0) is the initial angular position of the planet measured from the planet-moon barycenter. 
The quantities R, r, I, (fib and (fi p are shown in Figure [T] The functional form of TOA pertyPm 
can be derived from R s , the vector between the system barycenter and the pulsar, using: 



o jo 



Rs • ndt dt — TO A per t^ p + TO A per t )Pmi 



(3) 



where c is the speed of light and n is a unit vector pointing along the line of sight, the only 
direction along which quantities can be measured. 

The governing equation for R s can be written as the sum of the zeroth order term, 
which describes TOA pert , p , and the tidal terms, which describe TOA pert , pm : 

d 2 R< 



G(M P + M v 



dt 2 



R 3 

G(M P + M, 



R 



R 



GM n 



IR 



:(R 



GM„ 



|R + r. 



fR- 



(4) 



where the tidal terms have been collected into square brackets and noting that R s = — (M p + 
Mm) f(M 8 + M p + M m )R, r p = -M m /(M P + M m )v and r m = M p /(M p + M m )r, and where 
the vectors R s , R, r p , r m and r are also shown in Figure [T] Using the coordinate system in 
Figure [TJ it can be seen after some algebra that: 



R + r rn 
R + r t 



R cos (fib — 
R cos (fib + 



M p ' 

— r cos 4> v 

M p + M m Vp 

M m ' 

— ——r cos (fi v 

M v + M m Vp 



i + 



i + 



R sin (fib — 
R sin (fib + 



Mp . ' 

— t —— r sin (fi v 

M D + M m ^ 



M„ 



M p + M n 



-r sin i 



(5) 
(6) 
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where i and j are denned in Figure [TJ and i is the direction to the line-of-sight, projected 
onto the plane of the orbit. 

As the orbits are both circular and coplanar, we have that 0&(£) = n b t + 0^(0), and 
4> p {t) = n p t + P (O) where n& and n p are the constant mean motions of the two respective 
orbits. Substituting equations ^ and Q into the coefficients of the last two terms of 
equation Q, assuming r <C R and using the binomial expansion to order r 2 /R 2 gives: 



GM n 



IR 



GM n 



1 + 3 



M 2 

p 



M n 



[M m + M p ) R 
2 / 3 15 



cos 



(M m + M p ) 2 R 2 



"2 + 



cos 



-4>p) 



(7) 



GM„ 



IR 



GM P 



M m r 
' (M m + M p ) R 



COS((p b 



Ml 



15 



cos 



(M m + M p ) 2 R 2 V 2 ' 2 
Substituting equations (p]), d6]), (ul) and (J8| into Q gives, after simplification: 



G(M p + M m ) GM p M m r 

Iv + 



2 r 



i? 3 (M m + M p ) i? 4 

x (cos + sin b j) — 3 cos(0{, — (p 



3 
2 



15 



cos 



,) (cos0 p i + sin0 p j) 



From Figure [TJ it can be seen that: 



n = sin Ji + cos Ik. 



(9) 



(10) 



Substituting equations ^ and (10) into equation ([3]), gives 



TOA 



pert,pm 



sin I GM p M m r 2 



An 



2 COS (j) b + 



c(M m + M p ) i? 4 
15 

cos(30b — 20 p 



8(n fc - 2n p ) 



cos 



-20 p 



2n p ) 2 



'111 



The cos 06 term in equation (11 ) has the same frequency as the signal of a lone planet and it 
can be shown that it acts to increase the measured value of M p + M m derived from TOA pert , p 
by (3/4)(r 2 / R 2 )(M p M m / \M P + M m )). Consequently, this term can be neglected as it will 
be undetectable as a separate signal. Also, the edge of the stability region for a prograde 
satellite of the low-mass component of a high-mass binary can be approximated by 0.36r# 
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for the case of circular orbits, where r# = R [(M P )/(3M S )] 1//3 is the secondary's Hill radius 



(Holman & Wiegert 1999). When r is equal to this maximum stable radius n, p m 8n&. As 
nb <C n p is likely, we have that the denominators of the cos (0b — 20 p ) and cos(30b — 20 p ) 
terms will never approach zero. This, in addition to the assumption of zero eccentricities, 



means that resonance effects can be neglected. Consequently, equation (11) can be simplified 
by neglecting in the denominators, giving: 



TOA 



pert,pm 



- sin IGM p M m r 2 
c{M m + M p ) R 4 [32n 2 p 



15 



cos(0b - 20 p ) + -— cos(30b - 20 p ) 



32n2 



Writing n p in terms of r, using Kepler's law, gives: 



TOA 



pert,pm 



sin / 



M p M m 



R / r \ 5 



32 



cos 



+ ^ cos ( 3< ^ 



20 p ) 
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(13) 



From equation (13), we have that the size of the perturbation varies as [M m Mp/(M m + 



M p ) 2 ][r/R] 5 times the system crossing time, R/c. So, the best hope of a detectable signal 
occurs when the planet-moon pair are close to the parent pulsar, widely separated from each 
other, both quite massive, and very accurate timing data is available. 



Our result is consistent with a similar study done by Schneider & Cabrera (2006), who 



calculated the radial velocity perturbation on one component of a binary star system for the 
case in which the other component was an unresolved pair. Converting their radial velocity 
perturbation to a timing perturbation, setting the mass of the planet and moon equal to 



each other, as in the case investigated by Schneider & Cabrera (2006), and noting that r in 



this work is equivalent to their 2a a, our formula and that of Schneider & Cabrera (2006) 
agree. 



4. Is is Possible to Detect Moons of Planets Orbiting Millisecond Pulsars? 

To investigate whether or not it is possible to detect moons of pulsar planets, we simplify 



equation (13) by summing the amplitudes of the sinusoids, giving the maximum possible 
amplitude: 

(rrnA \ 9 sin J M m M p R / r \ 5 
max (TOA pert , pm ) = 16 {Mm + M p )^ Kr) " (M) 

Thus, for a given r/R, the maximum amplitude increases linearly with R. So, for planet- 
moon pairs that are far from each other and their parent pulsar, detection may be possible. 
For example, a 0.1AU Jupiter- Jupiter binary located 5.2AU from a host pulsar would produce 
a TOA pertjPm of amplitude 960ns, which compares well with the 130ns residuals obtained from 



one of the most stable millisecond pulsars, PSR J0437-4715 (van Straten et al. 2001). 
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Fig. 2. — The regions of parameter space containing detectable (shaded) and stable (cross- 
hatched) moons of the planet PSR B1620-26 b are shown as a function of planet-moon 
separation and moon mass. The total mass and the distance of the planet-moon pair from 



the parent pulsar were assumed to be 2.5 Jupiter masses and 23AU respectively (Sigurdsson 
et al. 2003), while it was assumed that sin / = 1. The mass of the host was set at 1.7 solar 



masses (the sum of the mass of the pulsar and its white dwarf companion). The 3a detection 
threshold was calculated assuming the ~ 40/is timing residuals given in Thorsett et al. (1999) 
are uncorrelated and that similar accuracy TOA measurements of PSR B 1620-26 continue 
to the present day. The stability region was estimated as 0.36 times the Hill sphere of the 
planet. 



To demonstrate this method, the expected maximum signals from a moon orbiting each 
of the four known pulsar planets was explored. It was found that in the case of PSR B1620- 
26 b, signals that are in principle detectable could confirm or rule out certain configurations 
of moon mass and orbital parameters (see Figure [2]). 

In the particular case of PSR B1620-26 b, the perturbation signal will not exactly match 



the signal shown in equation (13) due to the effect of its white dwarf companion. The effect 
of this companion was investigated as a side project and it was found that its effect was 
to introduce additional perturbations on top of the TOA pertp and TOA pert)Pm calculated. 
Consequently, the detection threshold represents an upper limit to the minimum detectable 
signal and thus the analysis is still valid. 
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Unfortunately, there are practical limits to the applicability of this method. They 
include discounting other systems that could produce similar signals, sensitivity limits due 
to intrinsic pulsar timing noise, and limits imposed by moon formation and stability. 

First, other systems that could produce similar signals need to be investigated. Possible 
processes include pulsar precession (see e.g., Akgiin, Link & Wasserman 2006), periodic 



variation in the ISM ( Scherer et al.fl997 ), gravitational waves ( Detweiler||1979 ), unmodelled 
interactions between planets (Laughlin & Chambers 2001) and other small planets. To help 



investigate the last two options, we plan on completing a more in-depth analysis of the 
perturbation signal of an extra-solar moon, including the effects of orbital inclination and 
orbital eccentricity. 

Second, the noise floor of the system needs to be examined. The suitability of pulsars 
for signal detection is limited by two main noise sources, phase jitter and red timing noise 



(for example, Cordes 1993). Phase jitter is error due to pulse-to-pulse variations and leads 
to statistically independent errors for each TOA measurement. Phase jitter decreases with 
increasing rotation rate (decreasing P) due to the resultant increase in the number of pulses 
sampled each integration. Red timing noise refers to noise for which neigbouring TOA 
residuals are correlated. Red timing noise has been historically modeled as a random walk 



Cordes 



1980 



Boynton et al. 


1972 


Groth 


1975 



Kopeikin 1997). Red noise is strongly dependent on P. This relationship 
can be understood from the theoretical standpoint that red noise is due to non-homogeneous 
angular momentum transport either between components within the pulsar (e.g. Jones|1990 ) 



or between it and its environment (e.g. Cheng 1987). To illustrate the effect of these two 
noise sources on TOA accuracy, an estimate of their combined residuals as a function of P 
and P is shown in Figure [3J For comparison, the values of P and P of every pulsar as of 
publication are also included. The ATNF Pulsar Catalogue^] ( [Manchester et aT 2005 ) was 
used to provide the pulsar data for this plot. 

Third, whether or not moons will be discovered depends on whether or not they exist 
in certain configurations, which depends on their formation history and orbital stability. 
Recent research suggests that there are physical mass limits for satellites of both gas giants 
( Canup fc Ward||2006 ) and terrestrial planets ( Wada fc Kokubo|[2b 06). Also, tidal and three- 
body effects can strongly affect the longevity of moons ( Barnes fc O , Brien||2002 Domingos, 



Winter & Yokohama 2006 Atobe & Ida 2007). 



Finally, while this method was investigated for the specific case of a pulsar host, this 
technique could also be applied to planets orbiting other clock-like hosts such as pulsating 



2 http: / /www. atnf.csiro.au/research/pulsar/psrcat/ 



-9- 



10 



-12 



-14 



-16 



10 



o 

ID 

« 10 

I 

60 -18 

a io 



10 



■20 



- 


■ " ■ ' y 


*4l»3K&«j^--. - 








r>^-" . 


















■ ■ 






/ 1 us| 10 us 


100 us 


1 ms 


10 ms 





10" 



10 



io" 1 

Log (P (sec)) 



10 



10 



Fig. 3. — Contour plot of predicted timing noise as a function of pulsar rotation period and 
period derivative. For comparison, all known pulsars as of publication are overplotted. This 



plot is based on Figure 9 from Cordes (1993). Consequently, the functions and assumptions 



used to generate the contours are the same as given in Cordes (1993), noting that that 
each TOA integration is 1000 seconds long. As the correlated timing noise measured for 
individual pulsars can vary from the predicted values by two orders of magnitude (see e.g., 



Arzoumanian et al. 1994), this plot is meant to demonstrate global pulsar properties, not 



predict individual pulsar noise characteristics. 



giant stars (Silvotti et al. 2007) and white dwarfs (Mullally, Winget & Kepler 2006). 



We are grateful to D. Yong for informing us of planet searches around pulsating stars. 
K. Lewis acknowledges the support of ANU/RSAA where the majority of this work was 
undertaken. 
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